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Abstract 

Based on the relation to random matrix theory, exact expressions for all microscopic spectral 
correlators of the Dirac operator can be computed from finite- volume partition functions. This is 
illustrated for the case of SU{Nc) gauge theories with Nc> 3 and Nf fermions in the fundamental 
representation. 
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Central to the understanding of chiral symmetry breaking in gauge theories is the spectral density p(A) 
of the Dirac operator. By the Banks-Casher relation p{0) = NT,/Tr , Tj = {ipip), this spectral density, 
when evaluated at the origin, is an order parameter for chiral symmetry breaking. As is well-known, 
this spontaneous symmetry breaking can, for the massless theory, only occur if one first considers the 
theory in a finite volume V, introduces a global source (mass) m for ij)^, takes the limit of 1/— >oo, 
and subsequently the limit m^O. One therefore needs to know the spectral density for small masses 
and large volumes, and then trace how the eigenvalues of the Dirac operator accumulate towards A = 
as the above limit is taken. 

A few years ago Leutwyler and Smilga greatly enhanced our understanding of this issue by deriving 
a series of exact spectral sum rules for the Dirac operator of gauge theories in a range 1/Aqcd ^ 
<^ I/ttt,^. Here Aqcd is a hadronic scale in QCD, V is the space-time volume, and is the pion 
mass. In doing this, Leutwyler and Smilga simultaneously derived exact analytical expressions for 
the finite- volume partition functions Z^^ of a variety of different theories in sectors of fixed topological 
charge z/. Based on a quite remarkable connection to random matrix theory ^ |5|, these exact 
statements about the Dirac operator spectrum are now understood in a completely new light. The 
first step in this series of events was the observation that the so-called microscopic spectral density of 
the theory with massless fermions, 

can reproduce all QCD spectral sum rules if computed in a large- random matrix model ensemble 
(the precise ensemble being uniquely dictated by the symmetries of the Dirac operator). One key 
to the understanding of this surprising fact is the by now proven universality, within random matrix 
theory, of all massless microscopic spectral correlators, and in particular the microscopic spectral 
density itself Q. An infinite sequence of spectral sum rules for the massless QCD Dirac operator can 
thus be reproduced by one single function, the microscopic spectral density defined above. While this 
is no proof that the microscopic spectral correlators are exact statements about QCD in the above 
limit, it is highly suggestive. It has also been put to test in a direct Monte Carlo evaluation of the 
microscopic spectral density (for SU{2), in the quenched approximation) 

Very recently, the above considerations have been extended to the case of massive fermions in QCD [Q] 
(see also Q). Here the appropriate limit is double-microscopic: both eigenvalues A and masses m must 
be considered on a scale of comparable magnification. Both for this case, and for QCD with massive 
fermions in (2 + 1) dimensions, the corresponding double-microscopic spectral correlators have been 
shown to be universal in the context of random matrix theory ||8|, |l^] . Moreover, the double-microscopic 
spectral densities satisfy non-trivial massive generalizations of the spectral sum rules |^, |l^, |Tll . 

One obvious question is the following. Given the fact that (double-) microscopic spectral correlators 
computed in a certain large-A random matrix theory can be used to "saturate" the spectral sum rules 
derived from finite- volume partition functions, are these spectral correlators themselves derivable from 
the finite-volume partition functions? There are already some results which would seem to indicate 
that this is not the case. For example, the finite-volume partition functions, and hence all spectral 
sum rules, for SU{2) gauge theory with 1 fermion and SU{Nc > 3) gauge theory with 1 fermion are 

^Here "QCD" is considered in the general sense of a (3 + l)-dimensional gauge theory with gauge group SU{Nc > 3) 
and Nf fermion species in the fundamental representation of this gauge group. (Generalizations to other theories are 
considered in ^ and Whenever we refer to the finite- volume partition function of QCD in what follows, it is 

considered in the above "mesoscopic" range of volumes. 
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identical. However, their microscopic spectral correlators are different Q. We shall nevertheless show 
in this paper that all information about the microscopic spectral correlators in a very precise sense is 
contained in the finite-volume partition functions. How the apparent counterexample above can be 
understood will become clear as we proceed. 

To be specific, let us consider the most interesting case, that of gauge group SU{Nc), Nc > 3, with 
Nf fermions in the fundamental representation. The coset of spontaneous symmetry breaking is here 
SU{Nf)L X SU{Nf)ji/SU{Nf), and the corresponding matrix models are those of the chiral unitary 
ensemble Q. We begin our considerations in the large- random matrix model language. Here 
partition function reads, for the sector of topological charge H 



,mNf 



where 



M 



dW YldetiM + imf) exp 
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f \ 
[ W j ■ 
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(2) 



(3) 



Here W is a rectangular complex matrix of size x (A^ + |i^|). In the large- limit the space-time 
volume V of QCD is identified with 2N. We integrate over the Haar measure of W. 

In terms of the eigenvalues A, of the hermitian matrix W'^W the partition function can be written 
(ignoring unimportant overall factors) 



Zu (mi, . . .,'mNf) 



N 



Nf 



^ poo y I J 

/ -^^ »=i V /=i 



mj) e-^^("») 



i-l 



detjjA^- 



(4) 



It is of course a symmetric function in mj. Moreover, one reads off that apart for the overall factor 
of Yl'm'jr, the partition function for Nf fermions in the sector of topological charge is equal to the 
partition function of the same Nf fermions plus additional massless fermions of zero mass, in the 
sector of zero topological charge. For simplicity, we set i^ = in what follows, so that there are no zero 
modes. 

A convenient expression for the spectral correlators of random hermitian matrices (the unitary ensem- 
ble) has recently been given in |12]. The appropriate generalization to the case of our chiral unitary 
ensemble with measure (^), where we are interested in the correlators of eigenvalues zi of M rather 
than those, Xi = zf, of W^W, is straightforward. The two-point correlator, the kernel, is 



{z,z ;mi, . . . ,miv ) 



-{V{z'^)+Viz'2)) 



N-1 



zz'll ^{z^ + m2)(z'2 + ^2) ^ P,{z^)P,(z'^) , (5) 



i=0 



where Pi{z'^) are the usual (orthonormal) polynomials associated with the above matrix model (see, 
e.g., ref. [^). They of course depend on all masses nif. The kernel can now be expressed as a 
normalized random matrix integral itself: 



{z,z ;mi, . . .,mNf) 



-f(V{z^)+V{z"')) 



'Uf Jiz"^ +mj){z'^ +171^) 



•^0 i=i \ f=i 



-W(A,)I 



detijA}-^ 
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The last integral is over (A^ — 1) eigenvalues only. However, in the large- A'^ limit we shall consider 
below, this distinction can be ignored. Thus, in the large-A'^ limit we have 



f 

Zq (mi, . . . ,mNf,iz,iz ) 

X JWJT, ^ ' (7) 

Zq (mi, . . . ,mNf) 

where the matrix model partition function in the numerator is evaluated for a theory corresponding 
to {Nf+2) fermions, of which two have imaginary mass. By means of the usual factorization property, 
all higher n-point spectral correlation functions are then also explicitly expressed in terms of the two 
matrix model partition functions Zq^^^ and Z^^^'^\ The spectral density corresponds to the two 
additional (imaginary) masses being equal: 

p^'^-t'\z;mi, . . . ,mNj>) = lim k'^'\z, z;mi, . . . ,mNf) ■ (8) 



We now consider the double-microscopic limit in which (" = zN2ttp{0) and pi = miN2TTp{0) are kept 
fixed as A^— >oo. In this limit the pre- factor exp[— (A^/2)(y(2;^) + ^(z'^))] becomes replaced by unity. 
Identifying T, = 2ttp{0), this is also the limit in which we can compare with the finite- volume partition 
function of QCD. 

What is the relation between the finite- volume QCD partition function Zq^^\pi, . . . ,pNf) and the 

matrix model partition Z^^\pi, . . . In the mesoscopic scaling region, which corresponds to 

the double-microscopic scaling regime of the matrix models, they should just be proportional, with 
a proportionality constant that is independent of the masses. This will then provide us with the 
sought-for relation between (double-) microscopic spectral correlators of the Dirac eigenvalues, and 
finite-volume partition functions. For the kernel, the master formula is 

/ ^0 {P'l,---,PNf> 

The double-microscopic spectral density is thus 

Ps iC; P'l, ■ ■ ■ , PNf) = C'ICIIKC +Pf) T]^) , (10) 

/ ^0 {pi,...,PNf) 

and the double-microscopic n-point correlation functions are given by 

p^^^\Cl,...,Cn;Pl,---,PNf) = det Kg^'\Ca,Cb; Pi, P-Nf) ■ (11) 

a,b 

The proportionality constant C is still undetermined, but there are several ways to fix it. One simple 
procedure is to use the matching between the microscopic spectral density p^^^\C; pi, . . . , pNj_.) as 

>oo with the macroscopic spectral density at the origin, or, in the conventional normalization 
l/vr. 

As a first example illustrating this, consider the microscopic spectral density p^g\C) of quenched QCD, 
which formally corresponds to Nf = 0. To find it, we need the finite- volume QCD partition function 
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for two massive fermions of degenerate (rescaled) masses i/i. This was evaluated analytically already 
in ref. p| and found to be, in their normalization, 



(12) 



where /n(x) is the nth modified Bessel function. The corresponding denominator in eq. (|^), is in 
this case just an irrelevant constant which can be set to unity. By means of the relation Iniix) = i'^Jn{x) 
between modified and ordinary Bessel functions of integer order, this gives 



pP{C) = c\c\ iMCf + JiiC? 



Requiring oo) = l/vr yields C = l/2, and hence 



pf(c) 
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ICI MCf + JiiCf 



(13) 



(14) 



the known result 0]. No explicit large- random matrix model computation is needed. But we can 
do more than this. The finite-volume QCD partition function for Nj fermions of arbitrary masses has 
recently been computed analytically by Jackson, §ener and Verbaarschot ||l^. The result is (dropping 
an irrelevant overall factor): 



A(^2 



where the Nj x Nj matrix A and A(/i^) are given by 



A^ 



i<j 



By repeated differentiation of the Bessel function relation 

Xl'n{x) = nln{x) + Xln+l{x) , 



(15) 



(16) 



(17) 



and by making use of the invariance properties of the determinant, we can replace the matrix A of 
(p!^) by a more convenient expression: 



A^ 



(18) 



For the numerator of eq. (|9|) we need the (A^/-l-2) x {Nf + 2) matrix A with two of the entries being 
imaginary. This means that 

Aij = {-dy-'jj-iiCi) for i = l,2, (19) 

and otherwise (for i>3) as in (|T8|). Also in the Vandermonde determinant of ([l^ ) it is convenient to 
separate out explicitly the terms arising from the imaginary entries. These terms are 

-lid -CD {uic! + f^})icl + f^}) 

while the remaining terms are identical to those coming from the denominator (and they hence cancel). 
For convenience we now pull out a factor of (-1) from every second column of the matrix A. This yields 
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an overall factor of (— l)!'^//^] where [x] denotes the integer part of x, and the matrix A is redefined 
accordingly (we call it B below). Putting these pieces together gives 



K^""^^ (C t-u u \-r (-l)[^//^]+^yClC2 det B 



where the {Nf + 2) x {Nf + 2) matrix B is defined by 



Bij = {dy-'Jj-iiCi) for i = l,2 

Bij = (-^,_2)^"'/j-i(/Xi^2) for 3 < i < TV/ + 2 , (21) 



and the NfxNf matrix A is as in (|l^). To find the corresponding double-microscopic spectral density, 
we make use of the Bessel relation 

^[x"J„(x)]=x"J„_i(x) (22) 

to get 

Ps (C,W,...,m,)-C 2n/(c2 + ^2) (23) 
where the {Nf + 2) x {Nf + 2) matrix B is defined by 

Bij = {Cy-^Ji-2{0 (24) 



and Bij = Bij for i^l. The general n-point correlators follow from eqs. ( |llD and (20). It finally remains 

to fix the constant C. We do this as before by the matching condition p~^^\c, ^oo, /xi, . . . , //at^) = l/vr. 
This gives 

C = (-l)[^//2] . (25) 

Substituting this into eqs. ( pOD and (23), the results agree with what has recently been obtained by 
an explicit computation in random matrix theory As mentioned earlier, the case 7^ can be 
extracted from the general formula for z^ = by setting v fermion masses equal to zero in a theory of 
Nf + v fermions. 

The case of the ordinary unitary ensemble, conjectured relevant for QCD with an even number of 
flavors in (2 + 1) dimensions |l^, can be treated in entirely the same fashion now that the finite- 
volume partition function is known explicitly [10|. There are also analogous relations between finite- 
volume partition functions and (double-) microscopic spectral correlators for the cases corresponding 
to orthogonal and symplectic random matrix ensembles. 



We end by a few comments on how these results can be intuitively understood. As we have seen, 
the (double-) microscopic correlators for the theory with Nf fermions can not be computed from the 
associated finite- volume partition function for Nf fermions alone. We need also to know the finite- 
volume partition function for two more flavors. These additional fermions, of imaginary mas^, act 
as sources that turn the partition function into a generating function. It is not surprising that the 
kernel, representing two degrees of freedom need two such additional fermions to probe the dynamical 
distribution of Dirac eigenvalues. To close full circle, we still need to understand how the master 
formula (^ can be derived directly from QCD alone. 

^There is no ambiguity associated with the analytic continuation to imaginary mass here. The finite-volume partition 
functions are explicitly given in terms of Bessel functions, defined in the whole complex plane. 
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